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ABSTRACT 

We investigate implications of quark-hadron duality for hybrid mesons 
in the large- iVc limit. A simple formalism is developed which implements 
duality for QCD two-point functions of currents of quark bilinears, with any 
number of gluons. We argue that the large- A'^c meson masses share a com- 
mon parameter, which is related to the QCD string tension. This parameter 
is fixed from correlators of conserved vector and axial-vector currents, and 
using lattice QCD determinations of the string tension. Our results predict 
towers of hybrid mesons which, within expected 1/ A'^c corrections, naturally 
accommodate the 1 experimental hybrid candidates. 
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1. Introduction 

The possible experimental observation of mesons with quantum numbers that cannot arise 
from a simple quark-antiquark state, has generated a great deal of interest [1,2]. In effect, 
search for these hybrid meson states is a primary mission of the proposed Hall D at 
Jefferson Laboratory. There are currently two experimental J^*^ — 1 ^ hybrid candidates, 
at ~ 1.4 GeV [1] and ~ 1.6 GeV [2]. One reason for excitement among theorists is the 
possibility that these states offer a direct experimental probe of the gluonic nature of QCD. 
Of course, a priori there is no reason why these meson states should not couple strongly to 
four-quark currents, rather than to currents involving gluons. In this respect, the large- A^c 
limit comes into its own as an extremely useful tool [3,4]. 

In the large- A'^c limit, mesons are narrow and couple to QCD currents with one quark, 
one antiquark and any number, Ug, of gluons [3,4]. We define a hybrid to be a meson 
whose space-time quantum numbers require Ug > 1. As with conventional mesons [3], the 
QCD correlators of currents which couple to hybrid mesons can be expressed as sums of 
meson tree graphs, which must be infinite in order to be dual to perturbative QCD at 
large momentum transfer [4]. By contrast, quenched lattice QCD calculations [5] suggest 
that the lightest hybrid has l""*" and a mass ~ 2 GeV. (A survey of results from hadronic 
models is available in Ref. 6.) The coupling of the hybrid mesons to currents with four- 
quark states is suppressed in the quenched lattice calculations. One might expect that 
the discrepancy between the lattice and experiment is due to this suppression [7,5]. The 
natural forum to test this hypothesis is the large- A^c limit, as the coupling of hybrid mesons 
to four-quark currents is subleading in the 1/A^c expansion [3,4]. 

In this paper we investigate duality in the large- A^c limit for QCD correlators with 
hybrid quantum numbers. We develop a simple formalism which implements duality for 
arbitrary two-point functions of currents of quark bilinears. The Wilson coefficients in the 
operator product expansion (OPE) of correlators with Ug = 1 also appear in the OPE 
of correlators with Ug = 0. We also argue that the large-A^c meson masses depend on 
a universal parameter. A, which is related to the QCD string tension. Recent work has 
investigated quark-meson duality in the large- A^c limit for conserved vector and axial-vector 
QCD currents [8,9] (see also Ref. 10). We use this phenomenology of ordinary mesons as 
input into our investigation of hybrids. 

In Section 2, we develop a simple formalism for investigating quark- meson duality in the 
large- ATc limit. In Section 3 we apply this formalism to conserved vector and axial- vector 
QCD currents and their well-known phenomenology. This section is a necessary review of 
work that has appeared elsewhere [9]. We discuss the issue of setting the scale at which 
to evaluate the QCD coupling constant in Section 4. We also introduce a universality 



3 



conjecture in this section. In Section 5 we investigate vector and axial-vector currents 
which couple to hybrid mesons. We discuss our results in Section 6. 

2. The General Statement of Duality 

The basic object we will consider in this paper is the two-point function, 

n (Q2) = i J d'^xe'i'^{0\T[3^ (x) J (0)]|0) , (1) 

where J{x) is a color-singlet QCD current of the form qT{ng) q, where T{ng) is an arbitrary 
Dirac-Lorentz structure coupled to Ug gluons, and = —q^. Here for illustration we have 
assumed a scalar current with respect to all quantum numbers. The general form of the 
spectral function of mass-dimension 2m in the large- A^c limit is 



-Im n (t) = 2 y f2 (n) m2"^ (n) 5 (t - (n)) , (2) 

n=0 

where F{n) is the decay constant of the nth meson, M(n) is its mass, and N regulates 
the number of mesons, which is strictly speaking, infinite. This spectral function expresses 
the statement that in the large- A^c limit, Il{Q'^) is given by an infinite sum of zero-width 
mesons. We can now construct Il{Q^) using dispersion theory. At large space-like momenta 
we obtain the general duality relation 

^ 7712 / N 



Q2 - (£^^))'^=2(^+™) 
Q2^^— - Q2i 



where C is an analytic function of ols which is computed in QCD perturbation theory 
and {0{Q'^,ii)) are Wilson coefficients of mass-dimension d = 2{£ + m). The a{N;iJ,)f^ 
are hadronic counterterms which render the hadronic sums over states independent of N, 
and is a renormalization scale defined in dimensional regularization with MS. The dots 
represent higher order perturbative contributions. (For a diagrammatic statement of this 
duality relation, see Fig. 1 and Fig. 2.) Eq. (3) implies 

F^{n)/M^{n) n'^. (4) 

n—^oo 

In this paper we will assume that the squared masses are linear in n. This is the case 
in two-dimensional QCD in the large- A/g limit [11]. And, too, linearity in n is expected 
if large- A^c QCD is in some sense equivalent to a string theory. There is a caveat to this 
assumption relating to chiral symmetry breaking, which will be addressed below. It follows 
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from Eq. (4) that F^{n) = F^ + 0{n-^) where F is a hadronic parameter and e is a positive 
number. We will keep the leading correction to F^. We choose the parameterization 



(n) = + A^n (5a) 

p2+2e 

Hence, a priori, each QCD two-point function is parameterized by four hadronic param- 
eters: M, A, F and E. It is straightforward to convince oneself that in order that the 
leading Wilson coefficient {d — 2m) not vanish, we must have e — m. On the hadronic 
side of Eq. (3) we then have 



2 (-1)"^ y ^ 2y-f—-^ + V a (iV; jj)^ Q^fc ^ 

- 2 ((-!)- g2. ^ + ^) (^^) - * (iV)) + ± a (iV; g2. ^q) 

where \1/ is Euler's function, defined as the logarithmic derivative of the Gamma function. 
The counterterms here are barred, indicating that they differ from those in Eq. (3). From 
Eq. (3) we immediately find 

C ii,) = i-ir^' ^ , (7a) 

Using Eq. (3) and Eq. (6), we can determine the counterterms, d{N] |U)fc, and match to the 
OPE separately for each m. Working to leading order in ag and keeping Wilson coefficients 
up to dimension six, we find 



7VA2 

a{N;iJ,)Q = C ill) log (8a) 

A* 

{o {^i)Y=' = c {^^) (m' - , (86) 

{O {fi)y=' = -\c ifi) (m^ - m2a2 + , (8c) 



{O (//))''=« = (//) (^m2 - iA2^ (m2 - a2) m2 , 
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for m = 0, and 



a (iV; = W (m) log ^ + ^ C [m^ - m2a2 + 1 , (9a) 
a(iV;/x)i = (^M2-iA2^ , (96) 

a{N;i^)^ ^ C{ix) log^, (9c) 

{o {^i)Y=' = Q c (/.) (m2 - a2) m2 + w (/.)^ (^m2 - Ia^^ , (9d) 

for m = 2. We will consider these two cases in this paper. 

One may wonder about the effects of higher-order contributions to F^{n). An 0{n~'^~'^) 
correction, parameterized by D^^^'^ / M'^'^'^'^ {n) , generates a logQ^ correction to the Wil- 
son coefficient with d = 2(1 + m). This correction is higher order in ctg than the order to 
which we are working. 

3. Conserved Vector and Axial- Vector Currents 

We consider QCD with two massless flavors. This theory has an SU{2) x SU{2) chiral sym- 
metry. In the large- A'^c limit, assuming conflnement, this symmetry breaks spontaneously 
to the SU{2)v isospin subgroup [12]. (Strictly speaking, large- A^c QCD has a U{2) x U{2) 
chiral symmetry. In this paper the additional Goldstone mode and its associated U{1) 
charge are ignored.) Consider the conserved QCD vector and axial- vector currents 



These currents fill out a six-dimensional chiral multiplet; they transform as (3, 1) ® (1, 3) 
with respect to SU{2) x SU{2). We are interested in the correlators 

2i J d'xe'^'^{0\Tmx)Vmm = - g^''q')ny{Q^) , (11a) 

2i I d'xe"i-{0\T[A^,ix)Aim\0) = SaMq"" - 9^"" q^)^A{Q'') • (12a) 

Since n^,A(<3^) are dimensionless, we are here interested in the case m = 0. 

Chiral symmetry constrains duality. Since chiral symmetry is not broken in perturba- 
tion theory, Ily{Q'^) = IIa{Q'^) to all orders in ctg. Moreover, because the pion couples to 
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Figure 1 : The statement of duality for the correlators of conserved vector and axial- 
vector currents. The crossed circles represent an insertion of a current. The thick line is 
a meson propagator and the filled-in square represents a hadronic countertcrm. The thin 
and wavy lines arc quark and gluon propagators, respectively. The small fillcd-in circles 
denote a perturbative interaction and the large filled-in circles denote a condensate. We 
have not included permutations. 



the axial- vector current, consistency with chiral symmetry requires that the pion, together 
with its chiral partners, be left out of the tower of states [8,9]. The minimal realistic ansatz 
for the spectral functions consistent with chiral symmetry is 



-ImUyit) = 2Fl5{t - Ml) + 2 V F^{n)5{t - M^{n)) , (13a) 
1 

-ImU^it) = 2Fl5{t) + 2Fl5{t - M^J + 2 J] Fl{n)5{t - Ml{n)). (136) 

Here we have extracted the lowest-lying vector and axial- vector mesons, p and ai, respec- 
tively. These states, together with the pion, satisfy spectral function sum rules. This 
is equivalent to the statement that tt, p and ai, together with an isoscalar Sq, fill out 
a reducible (1,3) © (3,1) ® (2,2) representation of the chiral group [9]. The parame- 
ters are then related by a single mixing angle, via F^^ = Fp sin Fa^ = Fjr cot and 
Mp = Ma^ cos 0. One can further show that each vector meson in the infinite sum must 
be paired with a degenerate axial-vector chiral partner so that pair-by-pair they fill out 
irreducible (1, 3) © (3, 1) representations of the chiral group. It follows that Fy = Fa = F, 
Ay = Aa = A and My — M^ = M [9] . It is intriguing that only the states outside of the 
tower exhibit mass splittings. 

Reading directly from Eq. (7a) and Eq. (8) we then find the system of equations 



(^v,A — 



127r2 ' 



{OYy=l = = 2F^ csc^ + Cy,A ( M' - ^A' 



etc 



{Ory:l = ^{G^.G,,) = -2F^Mlc8c' 



lfv,A ( 



6 



(14a) 
(146) 

(14c) 
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TH 


EXP 




TH 


EXP 




INPUT 


1465 ± 25 




44.1 ±0.6 


37.4 ± 1.2 


M„, 


1465 ± 25 


1640 ± 12 ± 30 




96 ± 2 


122 ± 23 


Mpii 


1903 ± 45 


1700 ± 20 


Mai 


1073 ± 7 


1230 ± 40 


Ma" 


1903 ± 45 


NO DATA 


F 


137 ±5 


NO DATA 




2258 ± 56 


2149 ± 17 




14.8 ± 0.02 


9 ± 2 [16] 


Ma'l' 


2258 ± 56 


NO DATA 


Ct-s {G ixvG ixv) 


0.063 ± 0.006 


0.048 ± 0.03 [17] 


Lio (10-3) 


-5.5 ±0.1 


-5.5 ±0.7 


A 


1172 ±43 MeV 





Table 1: Predictions compared to data using = 93 MeV, Mp = 770 MeV and Mpi = 

1465 ± 25 McV as input [14]. The theoretical errors are from Afp/. AU masses and decay 
constants are in MeV, Lio is dimensionless, is in degrees, the gluon condensate is in 
units of GeV* and the quark condensate (squared) is in units of 10"^ GeV^. 

{0)r = -f^c^siqqf = 2FlM^^csc^4>+\Cy,^ [m^ - (M^ - A^) ,(14d) 

{0)r = f = 2F2m,4^ + \Cy,. [m^ - (M^ - A^) m\ (14e) 

The Wilson coefficients and Cy^ have been computed from the diagrams in Fig. 1 by 
Ref. 13. We input F-j^ = 93 MeV, Mp = 770 MeV and the mass of the lowest-lying vector 
excited state, p'(1450): 1465 ± 25 MeV [14]. The solution to Eq. (14) with these inputs 
is given in Table 1. We have included predictions for the first few vector and axial- vector 
states in the towers of states. Of course the entire spectrum is predicted once M and 
A are determined. We have also included a prediction for the chiral perturbation theory 
parameter Lio, which is determined by the mixing angle (p [15]- The error quoted in the 
predicted quantities is from p'(1450). There are, of course, many additional sources of 
error, including 1/Nc corrections, which conservatively generate an error of ±3F7r (the 
delta-nucleon mass splitting). The predicted values of the quark and gluon condensates 
are somewhat large compared to the independent determinations of Ref. 16 and Ref. 17, 
respectively. 

4. Scale Setting and Universality 

A crucial issue of principle and practice is that of setting the scale at which the strong 
coupling constant and Wilson coefficients are evaluated. Intuitively one would expect that 
this scale, which can be taken to represent the onset of perturbative QCD, is a universal 
parameter, common to correlators with any given set of quantum numbers. Proving that 




® + (8«8) 

I n 
n 



Figure 2: The statement of duality for the vector and axial- vector hybrid correlators. 

The symbols are explained in Fig. 1. 

such a universal scale exists is highly non-trivial and we do not purport to do so here. 
Nevertheless, we will present a plausibility argument that such a universal parameter is 
present in large-A'^c QCD. Consider the basic relation between the hadronic description 
and the perturbative logarithm in large- A'^c QCD. From Eq. (3) and Eq. (5), neglecting the 
counterterms, we have 

^^Q2 + M} + nA^ qSoo " Aj^^'^iV^Af + •••• ^^^^ 

Here the subscript z is a collective meson quantum number and the dots represent terms 
that are suppressed by inverse powers of and Ni. The perturbative logarithm does 
not depend on the hadronic scale Mi. The natural scale at which to evaluate cts, which 
we interpret as the onset of perturbation theory, is ~ Aj. Naively, the Aj can be widely 
disparate. This renders choosing the physical scale at which to evaluate as channel depen- 
dent. Now consider the following argument. Regulating QCD with an ultraviolet cutoff 
Aqo implies 

A^ = A^iA? = iV^A^ = iVgAf = iV4A2 = . . . . (16) 

This is simply the statement that in perturbative QCD there is one scale. While physics 
is independent of the cutoff Aqo, seemingly physics depends on the Ni. Say we choose to 
regulate the sums with A^i = A^2- Eq. (16) then implies a new relation between physical 
quantities, Ai — A2, which is not present if instead we choose A^i ^ N2. This ambiguity is 
removed if 



Ai = A2 = A3 = A4 = . . . = A . (17) 

There is then a universal scale, A, which is shared by all large- A^c mesons. This common 
physical scale marks the onset of perturbation theory. In what follows we will treat Eq. (17) 
as a conjecture. Using the two-loop beta function in the large-AT^ limit in MS with 
Ams = 380 MeV, we find ^^(A) = 0.11, where we have taken A = 1172±43 MeV from the 
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TH (A,B) 


EXP 


TH (A,B) 


EXP 




INPUT 


1392 ±25 [1] 


INPUT 


1593 ±8 [2] 




1421,1543 




1787, 1884 






1990.2078 


loD.'-iiS [2] 


2394. 2407 


XO DATA 


Myll 


2445, 2588 


NO DATA 


2987,3104 


NO DATA 



Table 2: Predictions of 1 ^ hybrids using experimental candidates as input. On the left 
(right), we assume that the ~ 1.4 GeV [1] (~ 1.6 GeV [2]) candidate is the ground state. 
Theoretical numbers are from Eq. (22) using (A) predicted condensates from Table 1 and 
(B) quark and gluon condensates from Ref. 16 and Ref. 17, respectively (see Table 1). 

analysis of Section 3 (see Table 1). Wc will use this value of as in the next section. One 
may wonder whether this scale is sensible. If we assume an underlying string description, 
this scale relates to the string tension via A = yhia. From our analysis in the previous 
section, taking A from Table 1, we then infer ^/o" = 468 ± 12 MeV. This is in agreement 
with lattice QCD determinations of the string tension [18]. Lattice calculations of Mp/ y/a 
from the GFll [19] and UKQCD [20] collaborations yield = 440 ±15 ±35 MeV, where 
the first error is statistical and the second is systematic [18]. Of course the universality 
described here is precisely what one expects from Regge theory arguments, and, indeed, 
the observed Regge slopes imply = 400 — 450 MeV [18]. 

5. Vector and Axial- Vector Currents with Glue 

Consider the vector and axial-vector currents 

V« = ^.^A«G^,7'^iT«9, (18a) 
^«=^,gA«G>-75^T%, (186) 

where Qs is the QCD coupling constant, G^n is the gluon field-strength tensor and A" is 
an SU{Nc) generator. These currents transform as (3, 1) ® (1, 3) with respect to SU{2) x 
SU{2). Since these currents are not conserved, has nonvanishing matrix elements 
between 1"+ and 0''"''" states, and A^^ has nonvanishing matrix elements between l+~ and 
states. At leading order in the 1/A^c expansion, the 0''"''" and 1+" states also couple 
to currents with Ug = 0. Therefore we will have nothing to say about these states here. 
There are no obvious chiral constraints between the 1"+ and hybrid states. The 
nonvanishing matrix elements between the hybrid states and the vacuum are 
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Av Ap 

Figure 3: The left (right) panel gives Mv (Mp) as a function of (Ap) as given by 
Eq. (22). The shaded horizontal bars give the experimental 1 ^ events at ~ 1.4 GeV [1] 
and 1.6 GeV [2]. The solid curves use the predicted condensates from Table 1 and the 
dashed curves use the quark and gluon condensates from Ref. 16 and Ref. 17, respectively 
(see Table 1). The blocks over the solid curves are the universality predictions where A is 
taken from Table 1. They correspond to = 468 ± 12 MeV. The dotted curve in the 
right panel is evaluated with as = 0.13. 



{0\A''^\P') = 6'''FpMjpf,. (196) 

Here V denotes a 1 ^ hybrid and P denotes a hybrid. The 0"*""^ and 1"' states are 
denoted S and A, respectively. The relevant spectral functions are defined through 

8i J d^xe"i-{0\T[V!i(x)Vaom = 5ab{q^q' - g^'^ q^)Tly{Q^) + Sam^q^UsiQ^) ,{2Qa) 
8i I d^xe''i^{0\T[A^MAm\0) = Sabiq^q'' - 9^''q^)UA{Q^) + datq^q^^^^ 

Since n^AQ^) 

are dimension four, we are here interested in the case m = 2. The spectral 
functions in the large- A^c limit are 

1 

-ImUy{t) = 2y F|(n)M|(n)5(t-M|(n)), (21a) 

TV ^ — ^ 
1 

-ImUp{t) = Fj{n)Mj{n)S(t - MUn)). (216) 
Reading directly from Eq. (7a) and Eq. (9) we find the system of equations 



11 





TH 


EXP 




TH 


EXP 




1190 ±20 


1392 ±25 


Mp 


4180 ± 5 


NO DATA 


Av 


INPUT 




Ap 


INPUT 






1G70 ± 40 


loO.'-i ± S 


Mr- 


4.'-i4() ± 1() 


XO DATA 


Mv" 


2050 ± 50 


NO DATA 




4495 ± 26 


NO DATA 



Table 3: Predictions for 1 + and hybrids using universality. As input we use 
A,^ = Ap = A = 1172 ± 43 MeV (see Table 1), which corresponds to a string tension of 
= 468 ± 12 MeV. 



a.s ... oi 



C^ = -^, W^ = -^(G,.G,.), (22a) 

Cp = -^, >V, = ^(G,.G,.), (226) 

{0)t' = -^^«,(gg)2 = {hZy {Ml - ) m| + Wy^ (^M| - ^A^ ^ , (22c) 

{0)t' = -l^asiqqf = (m| - A|) M| + Wp^ (m| - ^A^ ^ . (22d) 

The Wilson coefficients and Cy p have been computed from tfie diagrams in Fig. 2 by 
Ref. 21. The masses of the lowest-lying hybrids, My p, require knowledge of Ay p. The 
purest predictions we can make use the 1 ^ experimental candidates as input to determine 
Ay. These results are given in Table 2. We give predictions using the condensates extracted 
in the analysis of Section 3, and from independent determinations [16,17] (See Table 1). 
The spread in values gives a sense of the theoretical error due to truncating the OPE. In 
Fig. 3 we plot My (Mp) as a function of Ay (Ap) as given by Eq. (22). We give curves 
for both sets of condensate parameters. The sensitivity to ag is insignificant for My. In 
Fig. 3 we show a curve for Mp evaluated with ag — 0.13. 

6. Discussion 

If we treat the 1"+ experimental candidate at ~ 1.4 GeV [1] as the lowest-lying state, 
then we predict the first excited state at ~ 2 GeV, as compared to ~ 1.6 GeV [2] (see 
Table 2). This discrepancy is larger than one would expect from a 1/Nc correction. If we 
treat the 1""*" experimental candidate at ~ 1.6 GeV [2] as the lowest-lying state, then we 
predict the first excited state at ~ 2.4 GeV. One must then understand the role of the 
~ 1.4 GeV [1] experimental candidate. We reiterate that in the large- A^c limit, mesons do 
not couple to QCD currents with four quarks. If such mixing is relevant, it corresponds 
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to a large subleading effect in the l/Nc expansion. Without assuming universahty, there 
is no prediction for the lowest- lying hybrid states. With Ay and Ap between 1 — 2 GeV, 
My ranges between 1.1 — 1.7 GeV, while Mp ranges between 4.1 — 4.7 GeV (See Fig. 3). 
Assuming universality, and taking Ky = Ap = A= 1172 ±43 MeV from the analysis of 
Section 3, we find a lowest-lying hybrid ~ 1.2 GeV, which underestimates the mass of the 
low-lying experimental candidate, and a first excited state ~ 1.7 GeV, which overestimates 
the mass of the heavier experimental candidate (see Table 3). However, the discrepancies 
are consistent with expected 1 / A'^c corrections. 
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